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1 introduction 



The quantum properties of the electromagnetic field apparently are infiu- 
enced by the presence of magnetodielectric media. For examples the spon- 
taneous decay constant of atomic systems [1]- [5], the Casimir effect [6] and 
the statistical properties of light are modified in the presence of a polarizable 
or magnetizable medium [T]-[9]. The quantization of electromagnetic field 
is considered usually in the presence of two types of media, the dispersive 
lossy magnetodielectric media and non-dispersive one. In dispersive magne- 
todielectric media there is a temporally nonlocal relationship between the 
polarization (magnetization ) field and the electric (magnetic) field [TU]-|13]. 
In this type of media, in order to inclusion the lossy effects, the medium is 
modeled by a collection of harmonic oscillators and both the electromagnetic 
field and the responsive medium are quantized, [Il]-[20]. These quantization 
approaches cover both isotropic and inhomogeneous anisotropic media. In 
these schemes the permittivity and permeability of the medium is calculated 
in terms of the classical parameters applied in the lagrangian or Hamiltonian 
of the total system. 

In the dispersionless dielectric media the relationship between the po- 
larization field and the electric field is temporally local. To quantize the 
electromagnetic field in the presence of such media, the medium is not quan- 
tized directly. The polarization effects of the medium is introduced, in the 
lagrangian density of the total system, only by the its linear and nonlinear 
susceptibility tensors [21]- [28]. In this case the Euler-Lagrange equations 
are the macroscopic Maxwell equations in the presence of the non-dispersive 
medium. This method has been developed for linear or non-linear isotropic 
dielectric media. In the dispersionless responsive media there is a class of lin- 
ear polarizable and magnetizable media that the electric polarization (mag- 
netic polarization) is related linearly to both the electric and magnetic fields. 
Generally, the electromagnetic properties of these media are described by four 
tensors of the second rank. This kind of media are known as bi-anisotropic 
magnetodielectric media. In this paper proposing a Lagrangian density a 
fully canonical quantization of electromagnetic field is represented in the 
presence of such media. It is worthy to point out that the form of Maxwell's 
equations in the curved space-time , written in the cartesian coordinates, in 
the absence of any responsive medium, are identical to the form of Maxwell's 
equations in fiat space-time in the presence of a bi-anisotropic magnetodi- 
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electric media P9| ISO] . Therefore the quantization method demonstrated in 
the present work is apphcable for the quantized electromagnetic field in the 
curved space-time in the absence of any responsive medium. The organiza- 
tion of this paper is as follows: 

In sec. 2 a Lagrangian for the electromagnetic field in the presence of a bi- 
anisotropic magnetodielectric medium is proposed and a classical treatment 
a bout the problem is achieved. In sec. 3 using the lagrangian introduced in 
section 2 a canonical quantization of the electromagnetic field in the presence 
of a bi-anisotropic inhomogeneous magnetodielectric medium is represented. 
In sec. 4 the quantization method is generalized in the presence of the atomic 
systems. In sec. 5 the spontaneous emission of an initially excited two-level 
atom located in a bi-anisotropic medium is investigated. Finally the paper 
is closed by a summary in sec. 6. 



2 Classical Lagrange equations 

In the non-dispersive linear magnetodielectric media the displacement field D 
and the magnetic field H are related to the electric field E and the magnetic 
induction field B as El] 



D(r,t) = e«(r)E(r,t) + e^^\r)B{r,t), 

H(r, t) = //«(r)E(r, t) + fi^^\r)B{r, t). (1) 

The electric (the magnetic) properties of the medium are described by the 
permittivities e^-^\r) , e^'^\v) ( the permeabilities fj,^^\r) , fi^'^\v)) , respectively. 
Generally, for bi-anisotropic media the permittivities and permeabilities may 
be appeared as four tensors of the second rank. There are some symmetry 
relations between the permittivities and permeabilities of a medium. The 
most important symmetry relations between the permittivities and the per- 
meabilities of a non- dispersive magnetodielectric medium are the Onsager's 
relations [21] 

A realization of the constitutive relations ([1]) and the symmetry relations 
(|2]) is the electromagnetic field in a curved space-time in the absence of any 
responsive medium. It is well known that Maxwell's equations in a curved 



3 



space-time in the absence of external sources, when the cartesian coordinates 
is used, are as the form SO] 



V-B = V-D = 

where the cartesian components of the fields D and H are defined by 

A = SijEj + (G X H)i Bi = fii.Hj + (G X E)^ (4) 

with 



^ij — l^ij — y 9 ^jJ — li2,3 

^00 

G. = -^ z = 1,2,3 (5) 

^00 

where g^j,^ is the space-time metric and g is its determinant. The defini- 
tions (jl]) for the fields D and H can be rewritten in the form of constitutive 
relations ([T]) with the tensors e'^^\ £^'^\ /x^^^ given by 

(1) _ , 9^^ 90a 90n 9l3m 



^ii V g ^iaB^mnj 

^ 9m 9my/-9 

(2) _ _ 9nj 90m 



(1) _ _ 90n 9im 



(2) 900 /^x 

^ij = ~^/^ 

where eijk is the 3-dimensional Levi-Civita symbol. It is easy to investigate 
that the tensors defined by ([6]) obey the Onsager's relation ([2]). Therefore, 
applying the cartesian coordinates, the form of Maxwell's equations in the 
curved space-time in the absence of a responsive medium is similar to the 
form of Maxwell's equations in the flat space-time in the presence of a bi- 
anisotropic magnetodielectric medium. The above discussions are sufficient 
motivation to investigate the quantization of the electromagnetic field in the 
presence of a bi-anisotropic magnetodielectric medium. 
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Using the constitutive relations ([T]) Maxwell's equations (E]) can be rewrit- 
ten as 

V • B = (7) 
VxE^-- (8) 

V ■ [£«(r)E(r, t) + £(2)(r)B(r, t)] = (9) 
V X [^«(r)E(r, t) + /i(2)(r)B(r, t)] = ^[e«(r)E(r, t) + e^'\v)Bir, t)] (10) 

The classical Maxwell's equations (|9]) and (ITOl) can be obtained as a conse- 
quence of the principle of Hamilton's least action using the Lagrangian 

L{t) = J £r£{r,t) = J £r{^e[f {r)E\r,t)E^ (r^t) (11) 
+Uf{r)E\r,t)B^{r,t) - \iifX^)B\v,t)E^ {v,t) - ^ fiff {r) B\r , t) {r , t)} . 



where the Einstein sum rule has been applied. In f lTTj) . E = — ^ — Vy^ and 
B = V X A. This means that the vector potential A and the scalar potential 
Lp constitute the degrees of freedom of the electromagnetic field. It is easy to 
show that the Lagrange equation for the variable f leads to the Gauss' law 
while, using the Onsager's relations the Lagrange equation for the 
vector potential A gives the Maxwell's equation f lTOj) . 

As it is clear in the present formalism the non-dispersive magneto dielec- 
tric medium is not quantized directly and the effect of the medium is appeared 
as the classical tensors and fj,^'^\ However in the quantization 

of the electromagnetic field in the presence of dispersive lossy media, the 
medium itself is quantized p^-[20]. 



2.1 Gauge fixing 

For a consistent canonical quantization of the electromagnetic field, we need 
the extra degrees of freedom to be eliminated from the above Lagrangian 
using some appropriate gauge conditions. Here we apply the gauge condition 

V-(£«(r)A(r,t))=0, (12) 
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where previously has been used by Glauber and et al [25]. Combination of 
this gauge with the Gauss' law leads to expression 



V?(r,t) = - J Sr' G(r,r') V ■ [^^'^r) V x A(r',t 



(13) 



for the scalar potential, where G(r, r') is the Green function satisfying the 
differential equation 

V ■ (£(^)(r) VG'(r, r')) = -5{y - r'). (14) 

Let us investigate the symmetry property 

G'(r,r') = G'(r',r) (15) 

for this Green function. To investigate this symmetry feature we consider 
a generalized Green's theorem for two arbitrary scalar functions V5i(q) and 
V?2(q) as 

^ (fq (ypi V ■ (e(^)V(/^2) - V ■ (£(')V(/^i)) 

= [ipi - V2 £^^^V(/5i) ■ h ds (16) 

which can be proved using the symmetry property of the tensor e^^\ that is 
efj^ = e^jl^ . Choosing v'i(q) = G(q, r) and V52(q) = G(q,r'), and regarding 
f lT^ and the boundary condition that the Green function should vanish on a 
surface located at infinity, the expected symmetry property ( 1T5|) is deduced. 

Let us expediently resolve a squared integrable vector field F(r) into two 
components, as F = F" + F^. The component F" is defined by 



Fllfr) 



J d\'VG{r,r')V ■ (£«(r') F(r')) (17) 



which is a conservative field. The other part F"*" is given by 

F^ = F-FII (18) 

which satisfies the relation V- (e:*^"'^^(r) F^(r)) = 0. we call F" and F"*" as e''^^- 
longitudinal and e^^^- transverse components of the vector field F, respec- 
tively. Conveniently, the tensor projection operators P^{r,r') and Pj||(r,r') 



are introduced as 



PMr,r') = 5,,5(r-r')-P, (r,rO. 



(19) 



In terms of these projection operators, one can rewrite the definitions ([T 
and ffTSl) as 



i = 1,2,3 



(20) 



J2 I rfVi^^(r,r')F,(r') 

.7 = 1 



i = 1,2,3. 



(21) 



Using the poisson's equation ( [T4|) and the symmetry property f|T5l) it is easy 
to show that the projection operators P-*" and P" satisfy the transversahty 
relations 



J = 1,2,3, 



i,m=l 



and 



E ^^t(r,r') = 



i = 1,2,3. 



(22) 



(23) 



The gauge condition ([T^ shows that — Vy? and — ^ are e*^^^- longitudinal 
and e*-^-*- transverse components of the the electric field E, respectively. This 
can be seen by operating P-*- and P" on the electric field, that is 



dt 

dxi 



J2 I d'r' Piiry)E,ir') 

J2 [ d'r' Pjl(r,r')i?,(r') 
i=i ^ 



(24) 
(25) 
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According to the Eq. fl23l) . one can obtain the transverse component of any 
squared integrable vector field F(r) as 

3 

F^ir) = E / ^^■(^') ^ = 1' 2' 3- (26) 

which clearly satisfies V ■ = 0. 



3 Canonical quantization 

Before beginning a canonical quantization scheme, the extra degrees of free- 
dom should be eliminated from the Lagrangian of the system applying the 
constraints f[T^ and f[T^ . This can be done by substituting the scalar po- 
tential Lp from f[T^ into the Lagrangian flTT]) and doing some integration by 
parts and using the antisymmetry relation e-j^ = — /^jP- Then we have 

m = I dh - 4)(r)^ [(e«(r))- s(^)(r) V x A(r,t 

-i/.;f(Vx A),(Vx A),l (27) 



j dh j rfV G'(r,r') V • \e'^^\r)W x A(r,t)] V ■ [£(')(r')V' x A(r',t)], 

Having the Lagrangian density of the system, denoted by £, as the inte- 
grand in ( 1271) . the z'th cartesian component of the canonical conjugate of the 
dynamical variable A can easily be calculated in a standard way as 

n.(r,.) = If (28) 



£"'(r) ^'^'g|'''' -^'"(r) [(£'"(r))"' ^'^'(rjV x A(r,i)],^ 



As it is seen, the canonical conjugate variable 11 is a purely transverse vector 
field , that is V ■ 11 = , while according to the gauge condition (fT2|) the 
vector potential A is a purely e'^^^-transverse vector field. Now the canonical 
quantization can be accomplished in a standard fashion by imposing the fol- 
lowing equal time commutation relations between the cartesian components 
of the conjugate variables 11 and A as 

[A(r,t) n,(r',t)]=i^|(r,r'). (29) 



8 



Because of the transversality relations 0221) and fl23p . it is clear that the 
commutation relations (12^ is compatible to the conditions V ■ 11 = and 
V • (e''^^(r)A) = 0. The time evolution of any dynamical operator related to 
the electromagnetic field, in the Heisenberg picture, can be obtained by using 
of the Hamiltonian of the system. The Hamiltonian of the electromagnetic 
field can be written in terms of A and 11 as ( see the Appendix A) 

H{t) = \j dV{[eW(r)],ri n,(r,t) n,(r,t) + 4^V) (V x A), (V x A),} 

+\ j [e(i)(r)],r.i |n,(r,t) [e^'\r)V x A]^. + [e^'\v)V x A]^n,(r,t)} 

+\ j d'r [e^'\r)]T^' [s^'\r)V x A], [e^''\r)V x A], (30) 

Any consistent quantization of the electromagnetic field should be able to 
give the Maxwell equations as the equations of motion of the electromagnetic 
field in the Heisenberg picture. It can be shown that the classical Maxwell's 
equation ffTOj) can be reobtained in Heisenberg picture if one uses the Hamil- 
tonian fl5U]) and applies the commutation relations (^U^. The Heisenberg 
equation for the vector potential gives (see the Appendix B) 

e«A = n + [(£«)-! ^ (3^) 



which coincides to the definition of the canonical momentum density in 
Also the Heisenberg equation for the conjugate variable 11 leads to (see the 
Appendix B) 



n 



V X (^^(^^(r) A(r,t)) -V X {fi^^\r)V x A(r,t)) + V x (/iW(r)V<^(r, t)) 

(32) 



where the Onsager's relations ([2]) has been used. Combination of the two 
Heisenberg equations (13T|) and (132!) yields 

= V X (;uW(r) A(r,t)) - V x {i2^^\r)V x A(r,t)) + V x (^«(r)V(/p(r, t)) 

(33) 

which, regarding the relation flA4p (in Appendix A), is precisely the Maxwell 
equation flTUl) . 
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3.1 diagonalization of the Hamiltonian of the system 

To obtain a diagonalized form of the Hamiltonian of the electromagnetic field, 
here for the sake of simplicity, we restrict ourselves to the special case that 
the tensors e^^-* and /i'-^-* are identically zero. Then, it is clear from Eq. f|T3|) 
that = and the Heisenberg equation fl55]) is reduced to 

- V X ^(2)(r)V X A(r,t) = e^'^\r)A{r,t). (34) 

Separating the time and space coordinates and writing A(r,t) = A{r)T{t), 
yields a differential equation for T as T = — oj^ T, where u"^ is a proper 
constant of separation, and the following equation for A(r) 

V X (r) V X A(r) = uh^^^ (r) A(r) , (35) 
In order to solve this equation, we consider a type of eigenvalue equation as 

V X /x(2)(r)V X F,(r) = ul e'^'\r) F,{r) (36) 

where is a squared integrable vector field that can be interpreted as an 
eigenvector field corresponding to an eigenvalue uf. The eigenvectors and 
the eigenvalues ul are obtained using the boundary conditions at infinity and 
the boundary conditions on the discontinuity surface of the medium. On the 
discontinuity surface of the medium, in the absence of external sorces, the 
tangential components of the fields E, H and the normal components of the 
fields D, B should be continuous. Using the symmetry property of the tensors 
e^^^ and given by ([2]), it can be shown that in the eigenvalue equation (l5^ . 
those eigenvectors which are correspond to the different eigenvalues, satisfy 
the orthogonality relation 

I d'r e[f{r)F*^ir) F,,,(r) ~ 5iu, - uy) (37) 

In order to show this. Let us write the complex conjugate of Eq. fl36|) for the 
label k' 

V X /.(2)(r)V X F^,(r) = u:*J e^'\v) F^,(r). (38) 

Muhiplying Eq.([3S]) by F^,(r) and Eq.([3S]) by Ffe(r) on the left, then sub- 
tracting and integrating the obtained results on the unbounded space and 
using the boundary conditions at infinity give us 

{^l - ^f) J d'r slf{r)F:,,{r) F,,(r) = 0. (39) 
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For k = k' the integral in the left hand of fl39p has a positive real value. To 
prove this, let C(r) be a complex tensor of the second rank defined by 

£(i)(r) = Ct(r) C(r) (40) 

where C^(r) is the Hermitian conjug ate of ^(r). Substitution e^^) from (00]) 
into the integral in (139|) . one can write 

I d'r e[f{v)F*,{r) F,,(r) = J d\ |C(r)F,(r)|2 > (41) 

provided that the tensor e^^^ and accordingly C(r) are assumed to be almost 
everywhere invertible. Equations ( 139|) and (HT!) clearly imply the reality of 
the eigenvalues uf. and the orthogonality relation (137|) . 

Although, applying ( |39|) and ( 14T]) . the orthonormality relation 

1 dhef^{v)Fl,{v) Fy,{v) = 8ky (42) 

can not be proved for those eigenvectors which are correspond to the degen- 
erate eigenvalues, in many cases it is possible to construct a complete set of 
the vector fields satisfying the eigenvalue equation f l5B]) and the orthonormal- 
ity relation ( 142|) . In fact in those cases that the tensor C(r), defined by (HOjl . 
is a hermitian tensor, that is = C, it can be shown that the differential 
operator 

K = [C(r)]-Wx/i(2)(r) Vx [C{v)]-^ (43) 
is a hermitian operator on the Hilbert space 

Vt = {g: ^ : j rf3rg*(r) ■ g(r) < oo} (44) 

where the inner product in f2 is given by 

Vf , gG^] (f|g) = j dVf*(r)-g(r) (45) 

Therefore there is a complete set of vector fields in f2 that satisfy the eigen- 
value equation 

[C(r)]-i V X /i(2)(r) V x [C{r)]-' h = ool h (46) 
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and the orthonormahty relation 



j (iVf*(r).ffc,(r) = 5fefc,. (47) 



It is clear from (l36l) and P6|) that the vector field Ffe(r) will satisfy Eg. (1361) 
iff C(r)Ffc(r) satisfies the eigenvalue equation (H6|l . Accordingly, at least, in 
those cases that the tensor C(r) is a hermitian tensor, one can construct a 
complete set of the squared integrable vector fields which satisfy the eigen- 
value equation fl36l) and the orthonormahty relation fH2|) . 

It should be noted that, as fl36|) shows, those eigenvector fields which are 
correspond to the nonzero eigenvalues satisfy the gauge condition 

V-[5(^)(r)Ffc(r)]=0. (48) 

Regarding to the completeness of the eigenvector fields one can expand the 
projection operator P-j-(r,r') ,given by ffH?]) . as 

^t(r,r') = EE ^^'(r') = EE ^^^w ^^(r )■ 

(49) 

where denote the summation over those eigenvector fields which are 

k 

correspond to the nonzero eigenvalues in fl36p . The condition (I48p shows that 
this expansion is compatible with the transversality relations fl22l) and (|23|1 . 
Also one can expand the vector potential A and the canonical momentum 
density 11 as 



A(r,t) = E \/2^ h(0)e-'^'^*F,(r) + 4(0)e^-'=* F*(r) 




hwh 




n(r,t) = ^ £(i)(r) E ^^^0) ^'^'^^ F^(^) " ^'^^^^ F'^'^^) 



(50) 

where and are the annihilation and creation operators of the bi-anisotropic 
non-dispersive magnetodielectric medium. Using the canonical commutation 



12 



relations fl29p and the expansions fl49p and flSOp . the commutation relations 
between and a|. is easily obtained as 

[afc(t), = Skk' (51) 

Regarding (HHj) . one can see that the gauge condition V ■ [e^^\r)A(r,t)] = 
and V ■ n(r,t) = are satisfied by the expansions (!50|) . Now inserting the 
expansions (HI]) in the Hamiltonian (I5U|) and applying Eqs.( 15B]) and fH21) . we 
reach to the diagonalized form for the Hamiltonian of the system (in the case 
e*^^-* = ji^^^ = 0) as the following 

i^F = ^ hwk a\ak (52) 

k 

where the normal ordering has been applied. 



3.2 Homogeneous bulk material 

For a homogeneous bulk material, that is when the tensors e*^^-* and //^^^ are 
independent of the position vector r, it is easy to show that the eigenvalues, 
ujI, in Eq. (!36|) are the roots of the following determinant 

det [A(q,//(2))_^2^{i)] =0, (53) 

where 

Aij(q,/i(^)) = -Eio^p ersj jJ'fr la Qs- (54) 

Let us denote the roots of the determinant in (1351) by ujp{q), where p labels 
the different roots of the determinant and q is an arbitrary three dimensional 
vector. Then, the normalized eigenvector field corresponds to the eigenvalue 
u!p{q) can be written as 

F(p, A, q, r) = _^^k^l^^_e^n-^ (55) 
v/Xt(p,A,q)£WX(p,A,q) ^ ' 

where X(p, A, q) is a three dimensional vector satisfying the algebraic eigen- 
value equation 

[A(q,/x(2))-a;J(q)e«]X(p,A,q)=0 (56) 

and A indicates the probable degeneracy of the eigenvalue ^^(q) which is 
known as the polarization of the photon. Here the label k in Eq.( l36i) is 
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represented by a triplet (p, A, q) and the orthonormality relation fl42p is ex- 
pressed as 

j £r e[f F*{p, A, q, r) F,(p', A', q', r) = 5,,,6xx'6iq - q') (57) 

The vector potential and the canonical momentum density 11 can be ex- 
panded in terms of the eigenvector fields (l55l) as 

^) = E E / ^""^^^ ^' ^' ""^ ^ ^'^^ 

n(r, t) = ^ J^JLJ^'i [at(p, A, q) e^-''('^)* F*(p, A, q, r) - H. 

Where denotes the summation over the nonzero roots in fl53|) . Finally, 

p 

the Hamiltonian of the electromagnetic field in the presence of the anisotropic 
homogeneous magnetodielectric medium can be written as 

P,A 

4 Generalization of the quantization in the 
presence of external charges 

The quantization method discussed in the previous sections can be general- 
ized in the presence of N point external charged particles. This is useful, 
particularly, when we are concerned with the spontaneous emission of the 
atomic systems within the magnetodielectric media. In the presence of A'^ 
charged particles with charges qi, q2, Qn and masses mi, m2, the 
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generalization of the Hamiltonian fl3Up is 

H{t) = \j dV{[e«(r)],ri n,(r,t) n,(r,t) + ^gV) (V x A), (V x A),} 
+\ J d'r [e«(r)],..i {n,(r,t) [e^'\r)V x A]^. + [e^'\r)V x A],^n,(r,t)} 
+1 j dh [e^'\v)]T^^ \s^^\v)V X A], [s^^\v)V x A], 

+E^^^^^^^ + E^^n-^'^) (60) 



i=l 



i=l 



where Xj and pj are the position and momentum operators of the ^'th particle 
that satisfy the usual commutation relations 



(61) 



where I is the unit matrix. In this case the commutation relations between 
the cartesian components of the conjugate variables A and 11 are the same 
as ( l29l) . In Hamiltonian ( |60l) . y(xj,t) is the coulomb potential produced by 
the other particles at the place of the i'th particle. This coulomb potential 
is related to the Green function introduced by f|T^ as 



(62) 



y(xi,t) = G'(xi,Xj) 



As the before section, in the special case that = 0, the conjugate 

variables A and 11 can be expanded in terms of the eigenvectors satisfying 
Eq.(l36]) as 



A(r,t) = 5^ 



h 




2uk 



afc(t)F,(r) + 4(t) F*(r) 



n(r,t) 



I e 



"WE 




4(t)F*(r)-a,(t) F,(r)l (63) 



From the commutation relations fl29p and the expansions 0491) . it is clear that 
the annihilation and creation operators a^, a|. satisfy the same commutation 
rules as fl5T]) . Now by inserting the expansions fl63|) into fl60l) and using the 
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eigenvalue equation (13^ . the generalized Hamiltonian (in the special 

case e^^-* = /i*^^^ = 0), is reduced to 



N 



[pi - g^A(x„f)]^ 
2m,; 



N 



+ y(x,t) (64) 



k i=l 



It should be pointed out that, in the Heisenberg picture, the time dependence 
of the operators ak{t),a\(t), appeared in the expansions fl63l) . is no longer 
sinusoidal as the before section. The time evolution of ak{t) and any other 
dynamical variable , in the Heisenberg picture, should be obtained by using 
the Hamiltonian (IMj) together with the commutation relations flHT]) and fl^T]) . 

5 Spontaneous emission of a two-level atom 

within an anisotropic magnetodielectric medium 

In this section, using the quantization method discussed in the before sec- 
tion, the decay rate of an initially excited two-level atom in the presence 
of an anisotropic non- dispersive magnetodielectric medium, is investigated. 
In the Hamiltonian flMj) . suppose that, except one of the particles, all the 
remainder have sufficiently large masses, so that they can be taken approxi- 
mately in fixed positions. Then, in the electric dipole approximation [32| 133], 
the Hamiltonian ( IMl) can be written as 



where e and m are the charge and mass of the particle that can move, re- 
spectively. 0(x) = e >^ G(x,Xj) is the coulomb potential at the place 




of the moving particle due to the other fixed particles. In (!65|) the position 
vector To is the center of the reign over which the moving particle is free to 
move. Such a collection of particles can constitute a one electron atom such 
that fq points to the center of the atom. In the Hamiltonian fl65|) the term 
^A^(ro) has been ignored , because in the electric dipole approximation 
this term dose not affect the decay rate of the atom [33j. To calculate the 
spontaneous emission of an atom, we restrict ourselves to the ideal model of 
a two-level atom. In this model the basis of the Hilbert space of the atom 



H = ^ hujk a\ak + h 0(x) p ■ A(ro) 




(65) 
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is contained two kets denoted by |1) and |2). The kets |1) and |2) are the 
eigenstates of the Hamiltonian of the atom corresponding to the eigenvalues 
El and E2, respectively. For the two-level atom the Hamiltonian f l65|) can 



now be rewritten as ^ 

H = Hp + Hat + H' 
Hp = ^ hujk a\ak 

k 

Hat = h 0(x) = hwo a'' a 

2m 

H' = HcjQfx d • A(ro) + H.C] (66) 

where a = |1)(2| , a"^ = \2){1\ are pauh operators, ojq = is the eigen- 

frequency of the atom and d is the electric dipole moment of the atom. Now 
by substituting the expansion f l63|) . for A(ro), into f l66|l the Hamiltonian in 
the rotating wave approximation reduces to |33j 

H = hujk a\ak + — iuJoe\^ \ [aa\ d • F^(ro) — H.C] (67) 

The simplest way to estimate the decay rate of an initially excited two- 
level atom is the Weisskopf-Wigner approach[33J. In this approach , in the 
schrodinger picture, the atom-field state at time t is taken as 

\m) = cm)\o) + j2 M,m)\k) (68) 

k 

where |0) is the vacuum state of the electromagnetic field and the the coef- 
ficients c{t) and Mk{t) satisfy the initial conditions c(0) = 1 , Mfc(O) = 0. 
Now applying the Hamiltonian (|67|1 and substituting \ipit)) from fl68|) in the 
schrodinger equation 

Him) = ^fi^^im) (69) 

for sufficiently large times, we have 

C{t) = e(-^'^«-'^'5-'-T)* (70) 

where 5u denotes the Lamb shift and 7 is the decay constant of the atom 
which are given by 

^-5i^5::^|d-F.(r.)r (71) 

k 
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7 = ^X.^'^ '^(^'^"^o) l^'^'^^'^o)!^ (72) 

k 

where P denotes the Cauchy principle value [25]. 

As an example let us assume that the center of the atom is located at 
the center of a very small spherical hole of radius R within the anisotropic 
magnetodielectric medium. Taking the origin of coordinates at the center of 
the hole, then the tensors e^^^ and fi^"^^ in spherical coordinates have the form 

£o r < R 

4V) = <! (73) 



f^i^)={ ■ (74) 
fiij r > R 

where Eq, /iq are the permittivity and permeability of the vacuum and Eij , fiij 
are two constant tensors which describe the electric and magnetic properties 
of the medium around the atom. Applying these tensors in the eigenvalue 
equation (l36l) . the following equation is deduced 

VxfiVxFk-ul eFk = Vx(/i- — ) e{R-r) VxFk-ul (e-leo) e{R-r) 

(75) 

where 6{R — r) is the unit step function and / is the unit tensor. If the 
radius R is assumed to be very smaller than the wavelength A = ^ then, 
the following approximations may be used for the i'the cartesian component 
of the right hand of fl75p 

[-ul {e - leo) e{R - r) Ffc(r)]^ ^ -ul {e,, - 6,,eo) e{R - r) F,,(0) 
V X (/i - — ) e(i?-r) V X Ffc 



Eij r > R 



j_ 

Mo 



r <R 



Mm - ^) -—m - r) F,,,„(0) + Q{R - r)F,,,„„(0) 



r 

(76) 

where Ffc,,„(0) and Ffc,,„„(0) denote |^(0) and (0), respectively. The 

solution of the equation fl75|) can now be written as the sum of two parts. 
One part is the solution of the homogeneous equation 

V X V X Ffc - a;2 £Ffc = (77) 
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that is 



Ffc(r) = F(p,A,q,r) 



X(p,A,q) 



(78) 



/Xt(p,A,q) £X(p,A,q) 

with X(p, A, q) given by (!53|) and (156!) . where the tensors e^^^ and /i*^^^ should 
now be replaced by e and /i , respectively. The second part of the solution 
(175]) is the response to the inhomogeneity term in the right hand of (175]) . The 
cartesian component of this part can be written as 



Fkiir) 



dh'G.,{v , r') J,(r') 



z = 1,2,3 



(79) 



where Jj is the sum of two approximated terms in ( 176]) and Gij{r,r') is the 
Green tensor related to Eg. (175]) satisfying 



(80) 



Using the technique of Fourier transform, the Green tensor G is easily ob- 
tained as 



G{y, r') = lim 



1 



»7-^o+ (27r)3 J 



1 



[A(p,p) - w2(q) e-iri\] 



(81) 



where / is the identity tensor and A(p,p) is defined by (!M]) with p*^^) = p. 
Using (176]) , (178]) and (179]) . the i'th cartesian component of the solution of 
Eg. (175]) can now be written as 



+ ^ 

oo 

+ 00 

+ 



^i(p, A,q) 



Xt(p,A,q) £X(p,A,q) 

+ 00 

Sr'G.jir, r') 9{R - r'){ej^ - 8^5^^) Ffc^(O) 



d r'Gij{r,r') ( ejap^mnsi^^p- 



6, 



e{R - r')Ffc,,„„(0) 



^6{R - r') Ffc,,„(0) 



(82) 



In ( 181]) and ( l82]) the eigenfreguency Uk has been replaced approximately 
by a;p(q), that is, the roots of the determinant [A(q, p) — w^(q) e\. Now 
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substituting (15T]) into ( 152]) . after integrating with respect to the variable r', 
it can be shown straightforwardly that the third term in (1821) is negligible for 
sufficiently small value of R and accordingly for the i'the cartesian component 
of the eigenfunctions Ffc(r), for <^ 1, one can write 



v/Xt(p,A,q)£X(p,A,q) 



n lim 

2n^ r;->0+ 



sm 



5. 



I3m ■ 



-^fcs,«7l(0) 

(83) 

Differentiation with respect to the cartesian coordinates x^, from the both 
sides of (1551) yields 



-ffcj,<57(r) 



v/Xt(p,A,q)£X(p,A,q) 



-qs q^)e 



— T lim 

27r^ r?->0+ 



+ 00 



d p -ps P-y 



Pks,an (0) 



Consistency condition at r = for the relations (183|) and (15^ gives 
^(1) p. ^n^ _ ^i(P>^>q) , ^(1) 



v/Xt(p,A,q) £X(p,A,q) 



isan ^ks,an 



(0) 



(85) 



and 



J2) 



Pks.an (0) 



Xi(p, A,q)(g5 g^) 



+ rgU^/^.«(o), (86) 



/Xt(p,A,q) £X(p,A,q) 

respectively, where the summation should be done over the repeated indices 
and the tensors V'^^^ , F^^^ , A*^-*^^ and A^^^^ are given by 

r> + 00 



r« = 5. +^lim 



d'p (^^) [A(p,/^)-c.2(q)e-^^/]--^ 



^7) 
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,(1) 



1 



o 9 1"^^ 



+ 00 



/Jo 



(2) .^^Mli^ 



+ 00 



+ 00 



+— lim 
27r2 T/-^o+ 



(90) 



From (l86|l one can obtain -Ffcs',Q,'n'(0) in terms of the cartesian component of 
F(0) as 



-Ffcs',o'n'(0) 



Xt(p,A,q)£X(p,A,q) 



+ 



where the tensor {^^'^^)s'a'n' is-y defined by 
Now, combination of fl85|) and fl9T|) yields 



(91) 



(92) 



Qjm-^fcm(O) 



VXt(p,A,q)£X(p,A,q) 



1 



where 



v/Xt(p,A,q) £X(p,A,q) 



n. - r^^) - A^^) rA(2)r^ r 



a£L (A^'^)s"i,,57 ^.(p>^'q) ^5 ?7 (93) 



(2) 



(94) 
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Finally, using the relation fl93p . one can write the m'th cartesian component 
of F{p, A, q, 0) as follows 



v/Xt(p,A,q)£X(p,A,(q) 



X 



qL"/) a, q) - gL7^ All ^.(P, a, q) gJ (95) 



Consequently, Applying the final result (|95|) in fl72|) . the decay rate of an 
initially excited two-level atom located at the center of a very small spher- 
ical cavity within an anisotropic homogeneous magnetodielectric medium is 
estimated as 

It can be shown, straightforwardly, in the case of an isotropic dielectric 
medium the result fl96l) coincides to the spontaneous emission of the atom 
computed by Glauber and et al in the limit R — )■ 0[25]. 



6 summary 

Proposing an appropriate Lagrangian density, a canonical quantization of 
the electromagnetic field in the presence of a non-dispersive bi-anisotropic 
inhomogeneous magnetodielectric medium was introduced. The quantiza- 
tion was achieved for both in the absence and in the presence of the atomic 
systems. The spontaneous emission of a two-level atom embedded in a bi- 
anisotropic homogeneous magnetodielectric medium was investigated. It was 
argued that, the form of Maxwell's equations in the curved space-time in the 
absence of a responsive medium , when the cartesian coordinates is used, 
are similar to the form of Maxwell's equations in the fiat space-time in the 
presence of a bi-anisotropic magnetodielectric medium. Therefore the quan- 
tization scheme discussed in this paper is applicable in the curved space-time 
in the absence of responsive media. This is useful particularly to investigate 
the effect of the curvature of the space-time on the quantum properties of 
the electromagnetic field. 
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Appendix A : Derivation of the Hamiltonian 

(SOD 

By definition, the Hamiltonian of the electromagnetic field is written as 

H{t) = j c/V[n(r, t) ■ A(r, t)] - L{t) (Al) 

where the Lagrangian L{t) is given by (1271) . If, using ( 128|) . the time derivative 
of the vector potential is obtained in terms of the dynamical variables 11 and 
A , then flAip reads to 

H{t) = \jd\ [{e^\'li, n, + /zSf (V X A), (V x A)," 
+i j dh |n ■ £(2)v X A]^ + £(2)V X A]^ ■ n} 

+1 J dh j dh' G{y, r') V ■ [£('Hr)V x A(r, t)] [V' • 5(')(r')V x A(r', t 

(A2) 

Because 11 is purely transverse, the second term is clearly equivalent to 

^ j dh {n . £(2) V X A] + £(2)^ X A] ■ n} (A3) 

According to the definition f|T7|) and the constraint f|T3|) one can write 

Vy^= £(2)V X A]" (A4) 

and therefore from f|T8l) the third term in flA2p is equal to 

(A5) 

which, after integrating by parts, is reduced to 

^ j d^r[e^%^^ [e^'^V X A]^ [e^^V x A], + ^Jd^ripV- [e^'V x A] 

(A6) 

Now applying Eq. (IT3|) the second term in flA6p cancels the last term in the 
Hamiltonian (IA2p and one can reach to the form given by fl5U]) . 
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Appendix B: Derivation of the Heisenberg equa- 
tion ([3ID 



According to the Hamiltonian (1301) the time evolution of the z'th cartesian 
component of the vector potential contains two parts. The first part is due 
to the commutator of Ai and the first term of the Hamiltonian (I5U1) . Using 
the symmetry feature of the tensor e^^^ and the commutation relations (129|1 . 
this part can be expressed as 

^ J d\' [e^'\v')]ll n„(r',t) [n,(r',t) , A(r,t)] 

= j <fr' [e^'\v')]ll P,i{r,r') U,{r',t) = [{e^'\r))-' Il{r,t)]^ (Bl) 

where the fact that the vector field (e*^^^)^^!! is purely 5^-transverse field has 
been used in the last step. The second part of the time evolution of Ai is 
caused by the commutator of the third and fourth terms of the Hamiltonian 
f pOj) and Ai, which using the symmetry property of the tensor e^^^ can be 
written as the following 

^ I dh' [e«(r')];^ {[H„(r',t) , A,(r,t)] (.(^)(r')V' x A(r',t))^} 

= jd'r' [e«(r)];^^i(r,r') (e^^^r') V x A(r', t))^ 

= VxA]^^ (B2) 

where the definition fl2Tl) has been applied in the last step. Now by adding 
fIBll) and flB2p the Heisenberg equation fl3Tl) is deduced. 

Appendix C: Derivation of the Heisenberg equa- 
tion (Bl) 

The Hamiltonian fl30l) shows that the time derivative of the i'th cartesian 
component of 11 contains three parts. Regarding the symmetry property 
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~ /^^a ^"^^ commutation relations fl29l) . the first part is 



2h 



d'r' /xS(r') \{V X A(r',t))„ (V x A{r',t)), , n.(r,t)' 



^ / /^S(r') (V X A(r',t))„ (V X A(r',t)), , n,(r,t) 



/ /.S](r') (V X A(r',t))^ e,,, a;P,i(r',r) 



(CI) 



Integrating by parts and applying ( l26i) . the final result in ( ICip is reduced to 
= - J d\' P,i(r',r) [V X V' x A(r',t) 



- [V X /i(2)(r) V X A(r',t 



(C2) 



Another part of Ilj is related to the commutator of the third and forth terms 
of the Hamiltonian (130|) and Ilj. Using the symmetry relations [^'-^^Jo^yg = 
[^'^^^]/3a ^"^^ commutation relations fl29l) . this part can be written as 



rfV [e^'\^')V, n„(r',t) [ {e^'\v')V x A(r',t))^ , n,(r,t) ] 

(C3) 

where the Onsager's relation e^^j = — /^^^ been used. Integration by 
parts and using the relations (|26|) , fl3T]) , ffTSl) and flA4l) the final result in fICSP 
is equivalent to 

= j d'r' P^^{v\v) [V'x/.«(r')(e«(r'))-'n(r',t)]^ 

= [Vx^W(r)(£«(r))-i n(r,t)]^ 

= [V X ;U«(r)A(r,t)], + [V x ^«(r)V(/.(r, t)], 
- [V X ii^^\T){e^^\r))-h^^\r) V x A(r,t)]. 

(C4) 
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Finally, the third part of time derivative of Ilj is caused by the commutator 
of the last term of the Hamiltonian (13U|1 and Ilj which can be computed as 

= y d?r' [e^'\v')V, {[ (^(^^(rOV' x A(r',t))„ , n.(r,t) ] {e^^Kv')V x A(r',t))^} 

= - / dh' [e^'\v')]ll e(2)(r') es„ d'Xdr',r) {e^'\r')V x A{r',t))^ (C5) 

Using the onsager's relation e^^j = — /U^^^, integrating by parts and then, 
applying the relation ( l26l) . this is equal to 

= j d^r' Pgi(r',r) [V x fi^^\r'){e^^\r'))-^ e^^\r') V x A{r',t)]^ 

= [V X //«(r)(5«(r))~i e^^\r) V x A(r,t)]. (C6) 

which cancels the last term in (IC4p . Therefore, adding the parts (IC2p ). flC4p 
and (IC6p gives the Heisenberg equation (l32l) . 
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